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Abstract 

We calculate 0(a s ) two- loop virtual corrections to the differential decay 
width dT(B — > X s £ + £~)/ds, where s is the invariant mass squared of the lep- 
ton pair, normalized to vn\. We also include those contributions from gluon 
bremsstrahlung which are needed to cancel infrared and collinear singular- 
ities present in the virtual corrections. Our calculation is restricted to the 
range 0.05 < s < 0.25 where the effects from resonances are small. The new 
contributions drastically reduce the renormalization scale dependence of ex- 
isting results for dT(B — > X s £ + £~)/ds. For the corresponding branching ratio 
(restricted to the above s-range) the renormalization scale uncertainty gets 
reduced from ~ ±13% to ~ ±6.5%. 
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I. INTRODUCTION 



After the observation of the penguin-induced decay B — > X s j and the corresponding 
exclusive channels such as B — > K*j ||, rare 5-decays have begun to play an important role 
in the phenomenology of particle physics. The measured decay rates are in good agreement 
with the standard model (SM) predictions, putting strong constraints on its various exten- 
sions. Another interesting decay mode in this context is the inclusive transition B — > X s £ + £~ 
{£ = e, /i). It has not been observed so far ||, but its detection is expected at the S-factories 
which are currently running. It is known that, unlike for B — > X s 7, large resonant contri- 
butions from cc intermediate states come into the game when considering B — > XJ + £r . 
When the invariant mass yfs of the lepton pair is close to the mass of a resonance, only 
model dependent predictions for these long distance contributions are available today. It 
is therefore unclear whether integrating the decay rate over these domains can reduce the 
theoretical uncertainty below ±20% 0. 

However, when restricting to regions of yfs below the resonances, the long distance effects 
are under control. In particular, all the available studies indicate that for the region 0.05 < 
s = s/ml < 0.25 these non-perturbative effects are below 10% ||- flQ| . Consequently, 



the differential decay rate for B — > X s l + l can be predicted precisely in this region using 
renormalization group improved perturbation theory. 

It is known that the next-to-leading logarithmic (NLL) result for the B — > X s £ + £~ 



decay rate suffers from a relatively large (±16%) matching scale (/j-w) dependence [|Tl~l , |T2" 
To reduce it, next-to-next-to leading (NNLL) corrections to the Wilson coefficients were 
calculated recently by Bobeth et al. |13| . This required a two-loop matching calculation of 
the effective theory to the full SM theory, followed by a renormalization group treatment 
of the Wilson coefficients, using up to three- loop anomalous dimensions [P^ , P^| . Including 
these NNLL corrections to the Wilson coefficients, the matching scale dependence could be 
removed to a large extent. 

However, this partially NNLL result suffers from a relatively large (~ ±13%) renormal- 
ization scale (fib) dependence (/i?, ~ 0(m,b)), as pointed out in ref. [|13|]. The aim of the 
current paper is to reduce this dependence by calculating NNLL corrections to the matrix 
elements of the effective Hamiltonian given in the next section. Our main contribution is the 
calculation of the 0(a s ) two- loop virtual corrections to the matrix elements of the operators 
0\ and O2, as well as the 0(a s ) one-loop corrections to O7-O10. Also those bremsstrahlung 
contributions are included which are needed to cancel infrared and collinear singularities in 
the virtual corrections. The new contributions reduce the renormalization scale dependence 
from ~ ±13% to ~ ±6.5%. 

The remainder of this letter is organized as follows. In Section II we review the theoretical 
framework. Our results for the virtual 0(a s ) corrections to the matrix elements of the 
operators Ox, O2, O7, 0$ and O9 we present in section III. Section IV is devoted to the 
bremsstrahlung contributions. The combined corrections (virtual and bremsstrahlung) to 
b — > s£ + £~ are given in section V. Finally, in section VI we analyze the invariant mass 
distribution of the lepton pair in the range 0.05 < s < 0.25. 
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II. THEORETICAL FRAMEWORK 



The most efficient tool for studying weak decays of B mesons is the effective Hamiltonian 
technique. For the specific decay channels b — > s£ + £~ (£ = /i, e), the effective Hamiltonian, 
derived from the standard model (SM) by integrating out the t-quark and the W^-boson, is 
of the form 
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where Oj are dimension six operators and Cj are the corresponding Wilson coefficients. The 



operators can be chosen as [13 
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where the subscripts L and R refer to left- and right- handed components of the fermion 
fields. We work in the approximation where the combination (V* s V ub ) of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements is neglected; in this case the CKM structure 
factorizes, as indicated in eq. (|l|). 

The factors 1/gl in the definition of the operators O7, O9 and Oio, as well as the factor 



l/g s present in Og have been chosen by Misiak [ nj in order to simplify the organization of 
the calculation: With these definitions, the one-loop anomalous dimensions (needed for a 
leading logarithmic (LL) calculation) of the operators are all proportional to g^, while 
two-loop anomalous dimensions (needed for a next-to-leading logarithmic (NLL) calculation) 
are proportional to g*, etc.. 

After this important remark we now outline the principal steps which lead to a LL, NLL, 
NNLL prediction for the decay amplitude for b — > s£ + £~: 



1. A matching calculation between the full SM theory and the effective theory has to be 
performed in order to determine the Wilson coefficients C, at the high scale fiw ~ 
mw, Tn t . At this scale, the coefficients can be worked out in fixed order perturbation 
theory, i.e. they can be expanded in g^: 

,4 
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At LL order, only cf^ is needed, at NLL order also C\ l \ etc.. While the coefficient 



(2) 

C) , which is needed for a NNLL analysis, is known for quite some ||15|| , C g 



(2) 



and C 



(2) 
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have been calculated only recently fl3[ (see also flql). 
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2. The renormalization group equation (RGE) has to be solved in order to get the Wilson 
coefficients at the low scale /i& ~ m^. For this RGE step the anomalous dimension 
matrix to the relevant order in g s is required, as described above. After these two 
steps one can decompose the Wilson coefficients C^/if,) into a LL, NLL and NNLL 
part according to 

CM = C?\» h ) + ^tff V & ) + Tf^C?^) + 0(g!) . (4) 



3. In order to get the decay amplitude, the matrix elements (s£ + £~\Oi(fib)\b) have to be 
calculated. At LL precision, only the operator Og contributes, as this operator is the 
only one which at the same time has a Wilson coefficient starting at lowest order and 
an explicit 1/ ' g 2 s factor in the definition. Hence, in the NLL precision QCD corrections 
(virtual and bremsstrahlung) to the matrix element of Og are needed. They have 



been calculated a few years ago |lT|jr2|] . At NLL precision, also the other operators 
start contributing, viz. 7 (fib) and Oi (fib) contribute at tree- level and the four-quark 
operators Oi,...,O e at one-loop level. Accordingly, QCD corrections to the latter 
matrix elements are needed for a NNLL prediction of the decay amplitude. 



As known for a long time fl? |, the formally leading term ~ (l/g 2 )Cg°\fib) to the ampli- 
tude for b — > s£ + £~ is smaller than the NLL term ~ (I / g 2 s )[g 2 s / (IQn 2 )] Cg\[i b ). We adapt 
our systematics to the numerical situation and treat the sum of these two terms as a NLL 
contribution. This is, admittedly some abuse of language, because the decay amplitude then 
starts out with a term which is called NLL. 

As pointed out in step 3), 0(a s ) QCD corrections to the matrix elements (s£ + £~\Oi([ib)\b) 
have to be calculated in order to obtain the NNLL prediction for the decay amplitude. In 
the present paper we systematically evaluate virtual corrections of order a s to the matrix 
elements of Oi, 2 , 7 , 0$, Og and O w . As the Wilson coefficients of the gluonic penguin 
operators O3, ...,Oq are much smaller than those of 0\ and O2, we neglect QCD correc- 
tions to their matrix elements. As discussed in more detail later, we also include those 
bremsstrahlung diagrams which are needed to cancel infrared and collinear singularities 
from the virtual contributions. The complete bremsstrahlung corrections, i.e. all the fi- 
nite parts, however, will be given elsewhere |[20|| . We anticipate that the QCD corrections 
calculated in the present letter substantially reduce the scale dependence of the NLL result. 



III. VIRTUAL CORRECTIONS TO THE OPERATORS O u 2 , O7, O s , AND Og. 

In this section we present our results for the virtual 0(a s ) corrections induced by the 
operators Oi, 2 , 7 , 8 , and Og. Using the naive dimensional regularization (NDR) scheme 
in d = 4 — 2e dimensions, both, ultraviolet and infrared singularities show up as l/e n -poles 
(n = 1,2). The ultraviolet singularities cancel after including the counterterms. Collinear 
singularities are regularized by retaining a finite strange quark mass m s . They are cancelled 
together with the infrared singularities at the level of decay width, taking the bremsstrahlung 
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process b — > s£ + £ (7 into account. Gauge invariance implies that the QCD-corrected matrix 
elements of the operators Oj can be written as 

(s£ + r \Oi\b) = F^(0 9 ) trec + PP(0 7 ) tree , (5) 
where (0 9 ) t ree and (0 7 ) tree are the tree-level matrix elements of O g and 7 , respectively. 



A. Virtual corrections to 0\ and O2 



The complete list of Feynman diagrams for the two-loop matrix elements of the operators 
0\ and O2 is shown in Fig. 1. Our calculation follows the line of |I1|[T9|j where the contribu- 



tions of O2 to the processes B — > X s ^ and B — > X s g have been evaluated. There, the results 
have been found as expansions in terms of powers and logarithms of the small parameter 
m 2 /m 2 . The central point of the procedure is to use Mellin-Barnes representations 



of certain denominators in the Feynman parameter integrals, as described in detail in refs. 
ITMCD!]. In the present case, however, we have an additional mass scale: q 2 , the invariant 



mass squared of the lepton pair. For values of q 2 satisfying < 1 and < 1, most of 

b c 

the diagrams allow a Taylor series expansion in q 2 and can be calculated in combination 
with a Mellin-Barnes representation. This method does not work for the diagram in Fig. 
|l]a) where the photon is emitted from the internal line. Instead, we applied a Mellin-Barnes 
representation twice. We will explain this procedure in detail in ref. [pop . The diagrams in 
Fig. |l|e) finally, we calculated using the heavy mass expansion technique f2l| . 

Using these methods, the unrenormalized form factors F^ 7 ' 9 ^ of 0\ and O2, as defined in 
eq. (H), are then obtained in the form 

^ (7 ' 9) = E <$2 * MS) H'k m W , (6) 

i,j,l,m 

where s = q 2 /m 2 and rh c = m c /mb. i,j,m are non- negative integers and I = + 

1/2, — i + 2/2, We keep the terms with i and I up to 3, after checking that higher order 

terms are small for 0.05 < s < 0.25, the range considered in this paper. 

The counterterm contributions are of various origin. There are counterterms due to quark 
field renormalization, renormalization of the strong coupling constant g s and renormalization 
of the charm- and bottom- quark masses. We stress that we use the pole mass definition 
for both, m c and m^. Additionally, we also have to take operator mixing into account. The 
corresponding counterterms to the matrix elements (CjOj) are of the form 

(ao^aE^^-)- with ( 7 ) 



Most of the coefficients aM 1 needed for our calculation are given in ref. [13|. As some are 



new, we list those for i = 1, 2 and j = 1, 2, 4, 7, 9, 11, 12 that are different from zero: 
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FIG. 1. Complete list of two-loop Feynman diagrams for b — ► 57* associated with the operators 
0\ and O2. The fermions (6, s and c quarks) are represented by solid lines; the curly lines represent 
gluons. The circle-crosses denote the possible locations for emission of a virtual photon. 






FIG. 2. Some Feynman diagrams for b — > S7* or 6 — > s^ + £~ associated with the operators O7, 
Og and O9. The circle-crosses denote the possible locations where the virtual photon is emitted, 
while the crosses mark the possible locations for gluon bremsstrahlung. See text. 
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On and O12, entering eq. (0), are evanescent operators, denned as 

On = (s L l^l^l^T a c L ) (c L7 ^7^ 7 ^T a 6 L ) - 16 0! 

Oia = (sl7„i7„27m3Cz) (c l7 w 7 M2 7 M3 &l) - 16 2 . 
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Before we give the result for the renormalized form factors, we remark that only diagram 
If) (and also its renormalized version) suffers from infrared and collinear singularities. As 
this diagram can easily be combined with diagram 2b) associated with the operator Og, we 
will take it into account in the next subsection, when discussing virtual corrections to Og. 



We decompose the renormalized matrix elements of Oj 



1, 2) as 



(s£ + £- \CVOi\b) = 0^ ( ~ ) ^ a; (0 9 }tree + ^ (Octree 
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with Og = Og and 7 = ^ 7 . The form factors F- 9 ^ and read (using L M = ln(fi/m b ), 
L a = hi(s))* 

r 9) ( 1424 16 . 64 r \ r 16 T T ( 16 32 A _ 2 \ r A 
Fl = U + 243 t7r + 27 Lc J " 243 ^ + 1 1215 " 135 J L "° 

+ ( 4 8 L^ + f 16 32 m A L s 3 256 L 2 + f (9) (12) 

+ 12835 " 315 mc J^ S + 176545 " 8505 mc ) L * S 243 L " + A ' (12) 



F 9 (9) = 
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+ (" 945 + 1^ ^) ^ ^ + (" 2^5 + ^) ^ + f V + ® ' <"> 



^ ) =^,+/? ) . (14) 

The analytic results for /-[ 9 \ /j- 7 - 1 , J^ 9 -*, and Z^ 7 - 1 (expanded up to s 3 and (ml) 3 ) are rather 
lengthy. The formulas become relatively short, however, if we give the charm quark mass 
dependence in numerical form (for the characteristic values of rh c =0.27, 0.29 and 0.31). We 
write the functions ffl as 

if = E k a\h3) ? Li (a = 1, 2; b = 7, 9; i = 0, 3; j = 0, 1). (15) 

The numerical values for the quantities k^(i,j) are given in Tab. I and II. 



B. Virtual corrections to the matrix elements of O7, 0$ and Og 

We first turn to the virtual corrections to the matrix element of the operator Og, con- 
sisting of the vertex correction shown in Fig. 2b) and of the quark self-energy contributions. 
The sum of these corrections is ultraviolet finite, but suffers from infrared and collinear 
singularities. The result can be written as 



( S £ + r \CgOg\b) = cf ( ) F 9 w (0 9 ) trec + F^(0 7 ) tTec 



a 



An 



(16) 



with Og = % Og and C { g ] = g (cf } + f^Cf } ) . The form factors F 9 {9) and F 9 {7) read (keeping 
terms up to order s 3 ) 



,(9) 



16 20 A 16 „ 2 116 ^ , 



(17) 



^(7) 



2 / 1 1 , 

-3 s ( 1+ 2 s+ r 



(18) 



6 



4 9) (o, 

*?>(!, 
M 9) (2, 

k?\3, 



M 7) (o, 

k?{2. 

*£>(3 : 
M ?) (3. 



-12.327 
-0.080505 

-33.015 
-0.041008 
-76.2 
-0.042685 

-197.81 
-0.039021 



= 0.27 

+ 0.13512 ? 

- 0.067181? 

- 0.42492 i 
+ 0.0078685 i 

- 1.5067? 
+ 0.015754 i 

- 4.6389 i 
+ 0.039384 i 



-11.973 
-0.081271 

-28.432 
-0.040243 

-57.114 
-0.035191 
-128.8 
-0.017587 



= 0.29 

+ 0.16371 ?' 

- 0.059691? 

- 0.25044? 
+ 0.016442?' 

- 0.86486? 
+ 0.027909? 

- 2.5243? 
+ 0.050639?' 



-11.65 
-0.080959 
-24.709 
-0.036585 
-43.588 
-0.021692 
-86.22 
0.013282 



= 0.31 

+ 0.18223?' 

- 0.051864? 

- 0.13474?' 
+ 0.024753 i 

- 0.4738? 
+ 0.036925 i 

- 1.3542? 
+ 0.052023 i 



-0.72461 

-0.26156 
-0.00017705 
0.023851 
0.020327 
0.42898 
0.031506 



0.093424? 







- 0.15008 i 
+ 0.02054? 

- 0.20313? 
+ 0.016606? 

- 0.099202? 
+ 0.00042591 i 



-0.68192 

-0.23935 
0.0027424 
-0.0018555 

0.022864 
0.28248 

0.029027 



0.074998 ?' 







- 0.12289 i 
+ 0.019676 i 

- 0.175? 

+ 0.011456i 

- 0.12783 i 

- 0.0082265? 



-0.63944 

-0.21829 
0.0053227 
-0.022511 

0.023615 
0.17118 

0.022653 



0.05885 ? 







- 0.10031? 
+ 0.018302? 

- 0.14836? 
+ 0.0059255? 

- 0.12861? 

- 0.0155 i 



TABLE I. Coefficients in the decomposition of f[ 9 ^ and f[ 7 ^ for three values of rh c (eq. ([15])). 

where the function /^j contains the infrared and collinear singularities. Its explicit form is 
(using r = (m s /mb) 2 ) 

/. nf = feL -(l + s + is 2 + -s 3 ) + 1 Wi- In r + - In r - - In 2 r . (19) 
J t 3 V 23 ; 3e 3 3 K J 

At this place, it is convenient to incorporate the renormalized diagram If), which has not 
been taken into account so far. It is easy to see that the two loops factorize into two one- 
loop contributions. The charm loop has the Lorentz structure of 0$ and can therefore be 
absorbed into an effective Wilson coefficient: Diagram If) is properly included by modifying 
Cf^ in eq. (Pf) as follows: 

Cf — Cf mod) = Cf + (cf + ^Cf ) H , (20) 
where the charm-loop function H reads (in expanded form) 



2835 



-1260 + 2520 ln(///m c ) + 252sm- 2 + 27rmJ 4 + 4s^ b . (21) 



In the context of virtual corrections also the 0(e)-part of this loop function is needed. We 
neglect it here since it will drop out in combination with gluon bremsstrahlung. Note that 
H = h{fh 2 c1 s) + 8/9 ln(/i/mb), with h defined in fT2]|T3[ . 
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TABLE II. Coefficients in the decomposition of and for three values of m c (eq. (|i~5|)). 



We now turn to the virtual corrections to the matrix element of the operator 7 , con- 
sisting of the vertex- (Fig. 2a) and self-energy corrections. The sum of these diagrams is 
ultraviolet singular. After renormalization, the result can be written as 
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with 7 = f- 7 and Cf ] = C 7 L> . The form factors Ff> and F^' ! read 
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Note that for these expressions the pole mass for rrib has to be used at lowest order. 

Finally, we give the result for the renormalized corrections to the matrix elements of Og- 
The corresponding diagrams are shown in Fig. 2c) and 2d). One obtains: 



(s£ + r \C 8 8 \b) = Cf > ( ) F 8 w (0 9 ) tree + F^"(0 7 ) tTee 
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with Cg '* = Cg~\ The form factors Fg 9 ^ and read (in expanded form) 

(9) 104 32 2 / 1184 40 2 \ /14212 32 2 \ 2 
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8 9 27 V 27 9 / V 135 3 / 
/193444 560 ,\ .. 16 r / „ ., 



IV. BREMSSTRAHLUNG CORRECTIONS 



We stress that in the present paper only those bremsstrahlung diagrams are taken into 
account which are needed to cancel the infrared and collinear singularities from the vir- 
tual corrections. All other bremsstrahlung contributions (which are finite), will be given 



elsewhere [20[] 



It is known [[jT],|T3] that the contribution to the inclusive decay width coming from the 
interference between the tree- level and the one- loop matrix elements of Og (Fig. 2b)) and 
from the corresponding bremsstrahlung corrections (Fig. 2f)), can be written in the form 



dr 9g fa em \ 2 G 2 F m\ le \V* s Vtb\ , ^ 2 , , ( /^(o)\ 2 ,.,\ fOQ , 

-w = UrJ — 48^ — (1 - s) (1 + 2s) ( 2 ^ ) -^j ' (28) 

where Cg ^ = ^ (Cg°^ + f^Cg 1 ^; the function u)g(s) = u)(s), which contains information on 



virtual and bremsstrahlung corrections, can be found in [|TT|,0. Replacing Cg by C, 



(0) 1 T ^(O.mod) 







(see eq. (P0"D) in eq. (28), diagram If) and the corresponding bremsstrahlung corrections 
are automatically included. 

Similarly, the contribution to the decay width from the interference between the tree- 
level and the one-loop matrix element of O7 (Fig. 2a), combined with the corresponding 
bremsstrahlung corrections shown in Fig. 2e), can be written as 
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The function lu 7 (s), which is new, reads 
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Finally, one observes that also the interference between the tree-level matrix element 
of 7 and the one-loop matrix element of 9 (and vice versa) lead to an infrared singular 
contribution to the decay width. We combined it with the corresponding bremsstrahlung 
terms coming from the interference of diagrams 2e) and 2f). The result reads 



dT 79 = /a em \ 2 G 2 F ml pole \V t * s V t -^ 

ds V AlT / 487T 3 V " " 7T 



79 _ { ®em\ ^ F m b,pole \ V ts V tb\ n ^ 2 i f ~\-Q /V^ (0)^(0) 



s) 2 12 ■ 2^u 79 (s)Re (cfcf) . (31) 



For the function 0*79(3), which also is new, we obtain 

co 79 (s) = -t In - ~ Li(i) - ^ rr 2 - I ln(i) ln(l - s) (32) 

12 + 75, M A , 2 s(3-2s) , ... 1 5-9s 

q — In 1 - S - - ^ s + — . 

9 s 9 (1 — s) 18 1 — s 



V. CORRECTIONS TO THE DECAY WIDTH FOR B -> 



In this section we combine the virtual corrections calculated in section III and the 
bremsstrahlung contributions discussed in section IV and study their influence on the decay 
width dY{b — > X s £ + £~)/ds. In the literature (see e.g. [0), this decay width is usually 
written as 



dT(b -> XJ+t 



ds 

A + 2s) 



a. 



An 



i ' n b,pole \ V t*s V tb 
487T 3 



s~icS 
°9 



+ 



°10 



+ 4 (1 + 2/s) pf J 2 + 12Re (C 7 eff C 9 eff *) 



(33) 



where the contributions calculated so far have been absorbed into the effective Wilson co- 
efficients C 7 S , C 9 S and . It turns out that also the new contributions calculated in the 
present paper can be absorbed into these coefficients. Following as closely as possible the 



'parametrization' given recently by Bobeth et al. ||13|| , we write 

a s (fi) 



Ceff 
9 



■eft 



c c 7 



Ceff 
10 



1 + 



7T 



u 9 (s) \{A 9 + T 9 h(m 2 c , s) + Ug h(l, s) + W 9 h(Q, s 



Air 



c[ 0) f} 9) 



^ 0) F 2 (9) + 4 0) F 8 (9) 



1 + 



7T 

Oi s (jj) 



Air 



(c[ 0) Ff 7) + Ci 0) Fi 7) +A^F, 



(7) 



7T 



u 9 {s) A 
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(34) 



where the expressions for h(m 2 c) s) and 0*9(5) are given in ]T3[. The quantities 0*7(5) and 
^12 8 ) on the other hand, have been calculated in the present paper. We take the numerical 
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a = 2 5 GeV 


u = 5 GeV 


u = 10 GeV 

Lai 1 V V.I V - v 


a s 


0.267 


0.215 


0.180 




-0.697 


-0.487 


-0.326 


cf 


1.046 


1.024 


1.011 


V ( ' (J 


(-0.360,0.031) 


(-0.321,0.019) 


(-0.287,0.008) 




-0.164 


-0.148 


-0.134 


(4 UJ ,4 iJ ) 


(4.241,-0.170) 


(4.129,0.013) 


(4.131,0.155) 




(0.115,0.278) 


(0.374,0.251) 


(0.576,0.231) 


«\^) 


(0.045,0.023) 


(0.032,0.016) 


(0.022,0.011) 




(0.044,0.016) 


(0.032,0.012) 


(0.022,0.009) 


<A (0} A (1) ) 

1^10 ' ^10 > 


(-4.372,0.135) 


(-4.372,0.135) 


(-4.372,0.135) 



TABLE III. Coefficients appearing in eq. (||) for u, = 2.5 GeV, /i = 5 GeV and fj, = 10 GeV. 
For a s (fj,) (in the MS scheme) we used the two-loop expression with 5 flavors and a s (mz) = 0.119. 
The entries correspond to the pole top quark mass nit = 174 GeV. The superscript (0) refers to 
lowest order quantities and while the superscript (1) denotes the correction terms of order a s . 



values for A 7 , Ag, Ai , T 9 , Ug, and Wg from flI5|, while C[°\ C% and = Cg 0,ef ^ are 



taken from [19 . For completeness we list them in Tab. III. 



When calculating the decay width (^), we retain only terms linear in a s (and thus in ujg 
and uj-j) in \Cf\ 2 and \Cf\ 2 . In the interference term Re (Cf Cf*) too, we keep only terms 
linear in a s . By construction, one has to make the replacements uog — > u-jg and ujj — > CU79 in 
this term. 

Our results include all the relevant virtual corrections and singular bremsstrahlung con- 
tributions. There exist additional bremsstrahlung terms coming e.g. from one-loop 0\ and 
O2 diagrams in which both, the virtual photon and the gluon are emitted from the charm 
quark line. These contributions do not induce additional renormalization scale dependence 
as they are ultraviolet finite. Using our experience from b — > 37 and b — > sg, these contri- 
butions are not expected to be large. 



VI. NUMERICAL RESULTS 



The decay width in eq. ([33|) has a large uncertainty due to the factor ml pole . Following 
common practice, we consider the ratio 

#quark (s) = — -— — , 35 

T(b — > X c ev) as 

in which the factor m\ pole drops out. The explicit expression for the semi-leptonic decay 
width Y{b — > X c eu e ) can be found e.g. in fl~3[ . 
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0.05 



0.1 



0.15 

s 



0.2 



0.25 



0.05 



0.1 



0.15 
s 



0.2 



0.25 



FIG. 3. (a) The three solid lines show the fi dependence of R qvL ark(s) when including the 
corrections to the matrix elements calculated in this paper; the dashed lines are obtained when 
switching off these corrections. We set rh c = 0.29. (b) -R qU ark(s) for rh c = 0.27 (dashed line), 
rh c = 0.29 (solid line) and rh c = 0.31 (dash-dotted line) and /i=5 GeV. See text. 



We now turn to the numerical results for -R quar k(s) for 0.05 < s < 0.25. In Fig. 3a we 
investigate the dependence of -R qU ark(s) on the renormalization scale \i. The solid lines are 
obtained by including the new NNLL contributions as explained in detail in section V. The 
three solid lines correspond to fi = 2.5 GeV (lower line), fi = 5 GeV (middle line) and \i = 10 
GeV (upper line). The three dashed lines (again fi = 2.5 GeV for the lower, /i = 5 GeV 
for the middle and fi = 10 GeV for the upper curve), on the other hand, show the results 
without the new NNLL corrections, i.e., they include the NLL results combined with the 
NNLL corrections to the matching conditions as obtained by Bobeth et al. ||13|| . From this 



figure we conclude that the renormalization scale dependence gets reduced by more than a 
factor of 2. Only for small values of s (s ~ 0.05), where the NLL /i-dependence is small 
already, the reduction factor is smaller. For the integrated quantity we obtain 

# q uark = / ds R qu ^ k (s) = (1.25 ± 0.08) x 10~ 5 , (36) 

J0.05 



where the error is obtained by varying \i between 2.5 GeV and 10 GeV. Before our correc- 
tions, the result was -R qU ark = (1-36 ± 0.18) x 10~ 5 |0|. In other words, the renormalization 
scale dependence got reduced from ~ ±13% to ~ ±6.5%. 

Among the errors on -R qU ark(s) which are due to the uncertainties in the input parameters, 
the one induced by m c = m c /mb is known to be the largest. We therefore show in Fig. 3b 
the dependence of -R qU ark(S) on m c . Comparing Fig. 3a with Fig. 3b, we find that the 
uncertainty due to m c is somewhat larger than the left-over /i-dependence at the NNLL 
level. For the integrated quantity i? quar k we find an uncertainty of ±7.6% due to rh c . 

To conclude: We have calculated virtual corrections of order a s to the matrix elements 
of Oi, O2, O7, 0$, O9 and 0\q. We also took into account those bremsstrahlung correc- 
tions which cancel the infrared and collinear singularities in the virtual corrections. The 
renormalization scale dependence of -R qua rk(s) gets reduced by more than a factor of 2. The 
calculation of the remaining bremsstrahlung contributions (which are expected to be rather 
small) and a more detailed numerical analysis are in progress | |20| . 
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